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Abstract. We present a quantum Monte Carlo study of the one-body density 
matrix (OBDM) and the momentum distribution of one-dimensional dipolar bosons, 
with dipole moments polarized perpendicular to the direction of confinement. We 
observe that the long-range nature of the dipole interaction has dramatic effects 
on the off-diagonal correlations: although the dipoles never crystallize, the system 
goes from a quasi-condensate regime at low interactions to a regime in which quasi- 
condensation is discarded, in favor of quasi-solidity. For all strengths of the dipolar 
interaction, the OBDM shows an oscillatory behavior coexisting with an overall 
algebraic decay; and the momentum distribution shows sharp kinks at the wavevectors 
of the oscillations, Q = ±2nn (where n is the atom density), beyond which it is 
strongly suppressed. This momentum filtering effect introduces a characteristic scale 
in the momentum distribution, which can be arbitrarily squeezed by lowering the atom 
density. This shows that one-dimensional dipolar Bose gases, realized e.g. by trapped 
dipolar molecules, show strong signatures of the dipolar interaction in time-of-flight 
measurements. 



1. Introduction 

Trapped quantum degenerate gases are currently offering the possibility of a thorough 
investigation of fundamental models of correlated quantum liquids, spanning a wide 
range of regimes from weak to strong interactions [1J. A particularly appealing aspect 
of cold atoms is indeed the tunability of the interparticle interactions. So far most 
of the experiments have been performed in a regime in which the interaction can 
be faithfully described as a contact s-wave interaction; tuning an applied magnetic 
field close to a Feshbach resonance allows to control continuously the magnitude and 
sign of the associated s-wave scattering length J2J. The use of Feshbach resonances 
allows in particular to suppress almost completely the contact interaction, leading to 
the observation of the effects of weaker magnetic dipole-dipole interactions for atoms 
which have a net magnetic moment, such as 52 Cr [31 EJ E]. On the other hand, the 
recent realization of a stable degenerate gas of heteronuclear molecules [6] appears to 



Off-diagonal correlations in a one- dimensional gas of dipolar bosons 



2 



pave the way to the exploration of the physics of quantum fluids with strong (and 
variable) electric dipole-dipole interactions; in fact the application of an electric field on 
the heteronuclear molecules, or the dressing of molecules with a microwave field, can 
induce a tunable electric dipole. 

On the theoretical side, the presence of the dipolar interaction has been shown to 
lead to many intriguing effects when the system is confined to two spatial dimensions, 
such as roton-like excitations in the liquid regime [7j, spontaneous crystallization [8], 
and even supersolidity in presence of an optical lattice [7] . In the case of dipolar bosons 
confined to one spatial dimension, in the absence of an optical lattice Mermin-Wagner- 
Hohenberg theorem forbids crystallization even at T = due to unboundedly strong 
quantum fluctuations. Hence the system can be described as a Luttinger liquid for all 
strengths of the dipolar interaction QUI [HI H2]- While it can develop arbitrarily 
strong diagonal correlations, resulting in a so-called super- Tonks behavior p^3l fl^] . such 
correlations are always decaying as a power of the distance. 

All the existing studies of the one-dimensional dipolar Bose gas have mainly focused 
on diagonal correlations and on the static structure factor [9], [TUl [121 ESI EH]- Yet 
the most natural observables in cold atom experiments are off-diagonal correlations, 
namely the one-body density matrix (OBDM) and its Fourier transform, the momentum 
distribution, which is obtained from time-of-flight measurements. In this paper we focus 
our attention on these quantities, which we compute numerically exactly by quantum 
Monte Carlo simulations. In particular, we show that off-diagonal correlations exhibit 
very strong signatures of dipolar physics. In fact the OBDM has a characteristic 
oscillatory behavior at the ordering vector Q of the Wigner crystal obtained in the 
classical limit. Such a behavior translates into a strong feature in the momentum 
distribution, namely a sharp kink at the Q vector, and a strong suppression of the 
momentum population for larger wavevectors. This effect of momentum filtering in 
the range [— Q,Q] is completely controlled by the density of the system, so that the 
momentum distribution can be arbitrarily squeezed in momentum space, although the 
system never attains condensation at finite density. A study of the two-body problem 
suggests that momentum filtering is a special feature of long-range interactions already 
at the few-body level, although its manifestation in the iV-body problem is much more 
pronounced. 

The paper is structured as follows. Section [2] describes the system Hamiltonian 
and the numerical method; Section [3] presents a study of off-diagonal correlations in 
the two-body problem with dipolar interactions, contrasting it with the case of contact 
interactions; Section [4] focuses on the OBDM as compared to predictions from Luttinger 
liquid theory; and finally Section [5] discusses the results concerning the momentum 
distribution. 
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2. Model and methods 

We investigate a system of N bosons confined to a one dimensional box of length L with 
periodic boundary conditions. The system Hamiltonian reads 

h 2 - d 2 

« = -^Ei+E^-*,). a) 

i=l 1 i<j 

In the following, otherwise specified, we will focus on the dipolar potential 

V ^ X - X ' )= \^f (2) 
Introducing the particle density n = N/L and the rescaled variables x~i = Xi/r (where 
r = mVd/ti 2 is an effective length associated with the potential), the Hamiltonian with 
dipolar interactions can be rewritten in the convenient form 



H(Ryd) = (nr f 



d 2 , 1 

9 + (nro) 



y~j /y»2 * ' ^^^^ I ™ ™ _ I 3 

^2 



(3) 



Here the Hamiltonian is expressed in effective Rydberg units, Ryd = nrjimr^f) 
Hence nro gives the dimensionless strength of the dipolar potential, and it contains 
a fundamental property of scale invariance, typical of power-law decaying potentials: 
rescaling the density n — > n/X and the effective length ro — > Aro, the physics of the 
system is unchanged (apart from a rescaling of the box length, L — > XL, and, accordingly, 
of the positions, x — > Xx). In the following we will compare the behavior of the dipolar 
system with that of the Lieb-Liniger gas, characterized by the (^-potential 

V s (x - x') = g 5(x - x') . (4) 

with repulsive nature, g > 0. The associated scattering length is given by a = 
—2h 2 /(mg). The dimensionless strength of the potential can be measured by (n|a|) _1 . 

We investigate the ground-state properties of the ID gas of dipolar bosons by 
quantum Monte Carlo simulations, based on the continuous-space Worm Algorithm, 
which allows one to obtain the one-body density matrix straightforwardly [TBI IT7] . 
While the method is strictly valid for finite temperatures only, the zero temperature 
behavior can be recovered for low enough temperatures. In addition to Monte Carlo 
simulations, we also use standard exact diagonalization to study the problem of two 
interacting bosons, whose results are presented in the next section. 

3. Two-body problem 

The study of the off-diagonal correlations of the two body problem turns out to be quite 
insightful in the perspective of the iV-body problem, as some fundamental traits of the 
latter are already exhibited by the former. We numerically evaluate the ground-state 
wavefunction 0(r) for the relative coordinate of the two-body problem on a L = 21 ring 
(where I is an arbitrary length unit), and the associated momentum distribution 
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Figure 1. Left panel: Occupation of the lowest momentum states as a function of the 
interaction strength for the ground state of two bosons interacting via a ^-potential. 
We observe that the weight lost at k = is distributed over all momenta. Right panel: 
Momentum distribution for varying strength of the S potential. 

We begin our analysis of the two-body problem by investigating the Lieb-Liniger 
two-body problem, whose results are shown in Fig. [TJ The evolution of the momentum 
distribution with increasing strength of the potential term (decreasing n\a\) shows 
a depletion of the k = contribution, and an increase of all the finite-momentum 
components, compatible with a gradual squeezing of the two-body wavefunction 0(r) in 
real space. In particular for all strengths of the interaction, the momentum distribution 
shows a characteristic decay as l/k A for all finite momenta - a behavior which is 
consistent with that of the large- tail of n(k) in the iV-body problem |18j . This means 
that an increasing interaction strength is simply imposing a global rescaling factor to 
the momentum distribution at finite k. 

The behavior in the dipolar case is radically different, as shown in Fig [2j First of all, 
in the limit nr Q — > the system maintains a finite occupation of non-zero momentum 
states, due to the fact that, in this limit, dipolar bosons reproduce the physics of the 
Tonks-Girardeau gas of impenetrable particles [T0l[TI] (see Sec. [5]for further discussion). 
Moreover, upon increasing the potential strength nr$ we observe not only a suppression 
of the k = momentum component, but also at finite momenta \k\ > Q — 2nn, at least 
for moderate potential strength. Hence the increase of finite momentum components, 
associated with the squeezing of the two-body wavefunction in real space, is highly 
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Figure 2. Left panel: Occupation of the lowest momentum states for the ground 
state of two bosons interacting via a dipolar potential. We observe that, for increasing 
interactions, the population lost at momentum k = is mainly transferred to 
momentum k — 7m, while the population at higher momenta can be even suppressed. 
This behavior is to be contrasted with that of contact interactions, Fig. [T] Right panel: 
Square modulus of the two-body wavefunction. 



inhomogeneous, and it only takes place within a momentum window < k < Q, 
while the momentum components out of this window are strongly suppressed. This 
inhomogeneous redistribution of populations in momentum space, driven by increasing 
the interaction, results in a peculiar momentum filtering effect, which is present in an 
even more dramatic fashion in the iV-body problem, as it will be discussed in Sec. [5] 

4. One-body density matrix 

The one-body density matrix (OBDM) 

C(x-x') = (^(x)ip(x')) (6) 

for a system of N dipolar bosons in a box of length L = Nl with periodic boundary 
conditions, has been investigated via quantum Monte Carlo simulations, for interaction 
strengths ranging between nr$ = 0.2 and nr^ = 5. The number of particles utilized for 
the results shown here is N=30, but we performed calculations with N=15 and N=60 
as well, in order to gauge the importance of finite-size effects. We generally find that 
simulations at a temperature k^T < 0.1 nr Ryd do not show any significant thermal 
effect, and all results presented in the paper are referred to this temperature range. 
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Fig. [3] shows the OBDM for two representative values of the potential strength nr^. 
Strong quantum fluctuations in one dimension strongly suppress correlations, so that 
the OBDM displays an algebraic decay. In addition to this behavior, we observe an 
oscillation with period 1/n, revealing the strong role played by interactions. Indeed, 
as we will see later, interactions lead to arbitrarily strong diagonal correlations, which 
imply a crystalline structure of the particles at short range. This structure is reflected 
in the OBDM, which exhibits a modulation with an amplitude that decays spatially in 
the same way as the density correlations. 

According to the general Luttinger liquid theory, density-density correlations in a 
boson liquid display a dominant decaying behavior at large distances in the form 



where K is the Luttinger liquid exponent; Q = 2im is the ordering vector of the classical 
Wigner crystal; and d(x\L) = L\ sin(7nc/L)|/7r is the cord function. On the other hand, 
the OBDM is predicted to exhibit the following (asymptotic) decay [19J: 



C(r) 



d{r\Ly/( 2K ) 



cos(Qr) cos(mQr 



d(r\L) 2K m d(r\L) 2m2K 



Hence, besides the dominant decay of the type \r\~ 1 ^ 2K \ the OBDM may also show 
a modulation at wavevector Q with an amplitude decaying as the density- density 
correlations, |r| _2is: , and modulations at higher harmonics mQ with increasingly faster 
decay. We can successfully fit our numerical results to Eq. [8j and we are able to 
resolve harmonics up to m = 2 within our numerical precision. In particular, for the 
coefficient b\ of the lowest harmonic we obtain negative values ranging in the interval 
b\ ~ —3 r5 x 1CT 2 . The negativity of this coefficient has important consequences for 
the shape of the momentum distribution, as will be discussed in the next section. 

In addition, fits to Eq. [8] allow to extract the exponent K, which is plotted in 
Fig. [4] as a function of the interaction strength nr$. K has the fundamental property 
of governing the scaling of the peak in the momentum distribution and in the static 
structure factor with increasing system size. Introducing the momentum distribution 

n(jfe) = \ I dr C{r) e lkr (9) 



Lj 

and the static structure factor 

S(k) = y [ dr G(r) e ikr (10) 
L Jo 

one obtains that n(k = 0) ~ £1-1/(2^ w hii e S(k = Q) ~ L^ 2A '. According to 
the general definitions, the system exhibits quasi- condensation when the peak in the 
momentum distribution diverges with system size (K > 1/2); and it exhibits quasi- 
solidity when the peak in the structure factor diverges (K < 1/2). As shown in 
Fig. |4j increasing dipolar interactions allow to continuously vary the K exponent from 
K = 1 in the weakly interacting regime (characterized by fermionization [TQl E]), to 



Off-diagonal correlations in a one- dimensional gas of dipolar bosons 



7 



K — » in the strongly interacting regime. In particular we find that for nr$ ~ 1.3 
the Luttinger exponent crosses the value K = 1/2, corresponding to a transition from 
quasi-condensation to quasi-solidity. Hence dipolar interactions allow to continuously 
tune the role of correlations, driving the system from a regime of dominant off-diagonal 
correlations to a regime of dominant diagonal ones. Our estimate of the K exponent 
stemming from the OBDM is in good agreement with the semi-analytic formula given 
in Ref. [TU] , extracted from a fit to quantum Monte Carlo results for the ground state 
energy. 
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Figure 3. One-body density matrix for N = 30 dipolar bosons. The dashed lines are 
fits to the Luttinger liquid formula, Eq. pj, truncated to the m = 2 order. 



5. Momentum distribution 

The evolution of the momentum distribution, Eq. [9j upon changing the interaction 
strength nr$ is shown in Figs. [5] and [6j As discussed in the previous section, the OBDM 
exhibits a modulation term at wavevector Q with a negative prefactor. This means that 
the corresponding feature in the momentum distribution is to be expected as a dip at the 
same wavevector. Indeed, the momentum distribution shows a sharp kink at k = Q, but 
it also shows a strong suppression of the momentum population at all momenta \k\ > Q, 
observed for all investigated values of nr . Hence we observe that the dipolar interaction 
introduces a characteristic scale Q in the momentum distribution, and a momentum 
filtering effect outside of the interval [-Q, Q\. This is to be contrasted with the case of 
contact interactions: indeed the Lieb-Liniger model has a momentum distribution with 
power-law wings, n(k) ~ k~ A [IE], and hence it does not contain any special momentum 
scale. The momentum filtering effect has the characteristic feature that, upon increasing 
the interaction strength nro, the weight lost in the k = peak is redistributed almost 
exclusively over the [— Q, Q] interval, similarly to what we have observed in the case of 
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Figure 4. Luttinger liquid exponent K as a function of the interaction strength utq. 
The shaded region marks the transition from quasi-condensation to quasi-solidity. 



the two-body problem. As a result, the strongly peaked distribution at low nr evolves 
into an almost triangular distribution for stronger nr$. For extremely large values of nr$ 
one might expect that the weight lost at k = will eventually start redistributing over a 
broader interval than [-Q, Q\. Nonetheless the persistence of Luttinger liquid behavior 
in the system for arbitrary large values of nr Q leads us to conclude that the characteristic 
periodic modulation of the OBDM, responsible for the suppression of the momentum 
distribution around q = Q, will also persist to much larger values of nr^ than the ones 
considered here, so that the momentum distribution will maintain a strong signature of 
dipolar physics around that momentum value. 

Momentum filtering can be probed experimentally in two different setups. Working 
at constant density n, nr^ can be controlled by increasingly polarizing the dipoles and 
hence increasing r . This leads to the observations of Fig. [5j in which the momentum 
distribution is always defined on the same support and it changes shape from strongly 
peaked to triangular. A second type of experiment can be carried out at fixed dipole 
polarization, by changing the dipole density n, e.g. by varying the trapping potential 
which holds the system. Fig. [6] illustrates such a protocol, which allows to shrink 
arbitrarily the support [-Q, Q] of the momentum distribution by lowering the trapping 
frequency, and hence to concentrate all the atoms to an arbitrarily small volume in 
momentum space. Yet paradoxically the system never condenses; quite oppositely, as 
discussed in Refs. [10J and [H], it fermionizes in the limit nr — > 0, reproducing the 
behavior of a Tonks- Girardeau gas. 
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Figure 5. Momentum distribution as a function of the interaction strength nr for a 
system with fixed density n = l/l. The dashed lines mark the momenta Q = ±2nn. 
The right panel shows the distribution in logarithmic scale, to evidence the sharp 
suppression for momenta |fe| > Q. 




Figure 6. Momentum distribution as a function of the interaction strength nro for a 
system with fixed dipoles, namely a constant tq. Here the number of particles is fixed 
to N = 30, and the system size changes when changing nro, and it goes from L = 150/ 
for nrg = 0.2 to L = 61 for nro = 5. 



6. Experimental realization 

The momentum distribution is notoriously the most accessible quantity in trapped gas 
experiments, simply obtained by time-of-flight measurements [I]. More challenging is 
the realization of a quantum degenerate gas with dominantly dipolar interactions [7]. 
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Quantum degeneracy has been achieved in 52 Cr (3j HI [5], for which s-wave interactions 
can be suppressed by a Feshbach resonance in favor of residual magnetic dipolar 
interactions. In this system one has Vd = fiof^m/i^ 71 ) where \i m = (Bohr magneton) 
is the magnetic moment of 52 Cr, and consequently ro = 2.7 x 10~ 9 m. For characteristic 
condensate densities of 10 13 cm -3 we obtain linear densities n ~ 2 x 10 6 m _1 , and 
an effective potential strength nr$ « 0.5 x 10~ 2 , which does not allow to explore the 
transition from quasi-condensation to quasi-solidity described in this paper, but it might 
be already sufficient to observe effects of momentum filtering. Polar molecules, on the 
other hand, can potentially feature dipolar interactions with characteristic lengths r 
which are 3 orders of magnitude higher than those for 52 Cr [7], and they would hence 
allow to fully span the regime of interactions described in this paper. 

Another potentially challenging issue is the confinement of dipolar gases in one 
dimensional traps. Indeed confining the gas in an array of tubes created by a two- 
dimensional optical lattice leads to g^asz-one-dimensional systems, with sizable residual 
inter-tube coupling due to the long-range nature of dipolar interactions. A possible 
solution is the use of largely spaced tubes, obtained via two-dimensional optical lattices 
created by crossing lasers at an angle far smaller than tt. Another solution is trapping 
in single-mode tight atom waveguides created e.g. by atom chips [20] . 

7. Conclusions 

In this paper we have investigated the effect of dipolar interactions on the off-diagonal 
correlations of a one-dimensional Bose gas. We have shown that the long-range nature 
of dipolar interactions changes radically the momentum distribution with respect to the 
case of contact interactions, introducing a characteristic momentum scale beyond which 
the momentum distribution is strongly suppressed. The effect of momentum filtering is 
fully controlled by the gas density, and it allows to squeeze the momentum distribution 
to an arbitrarily small portion of momentum space. The momentum distribution 
can also reveal the interaction-induced transition from quasi-condensation to quasi- 
solidity, via the loss of scaling of the zero-momentum peak with system size. All the 
phenomena discussed in this paper can be observed via time-of-flight measurements on 
polar molecules trapped in largely spaced two-dimensional optical lattices, or in tight 
atom waveguides. 
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